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Abstract 

We compute the simplest hard thermal loops for a spatial 't Hooft loop in the deconfined phase of 
a S'U{N) gauge theory. We expand to quadratic order about a constant background field Aq = Q/g, 
where Q is a diagonal, color matrix and g is the gauge coupling constant. We analyze the problem in 
sufficient generality that the techniques developed can be appHed to compute transport properties 
in a "semi" -Quark Gluon Plasma. Notably, computations are done using the double line notation 
at finite N. The quark self-energy is a Q-dependent thermal mass squared, ~ g^T"^, where T is 
the temperature, times the same hard thermal loop as at Q = 0. The gluon self-energy involves 
two pieces: a Q-dependent Debye mass squared, ~ g^ T^, times the same hard thermal loop as for 
Q = 0, plus a new hard thermal loop, ^ g^T^, due to the color electric field generated by a spatial 
't Hooft loop. 

PACS numbers: 11.10. Wx, 12.38.Mh, 25.75.-q 



I. INTRODUCTION 



The spectacular success of the experiments at the Relativistic Heavy Ion Colhder (RHIC) 
at Brookhaven has invigorated the study of gauge theories at a nonzero temperature, T 
[1, 2]. Experiment clearly shows that an understanding of quantities in thermal equilibrium 
are not sufficient to understand the data, and that one also needs quantities near thermal 
equilibrium, especially transport coefficients. 

In this paper we compute for what appears to be an unrelated problem: the real time 
response functions for a spatial 't Hooft loop [3, 4]. A Wilson loop, ti V exp{ig A^dx^), 
represents the propagation of a test electric charge along the path C, and measures the 
response to magnetic flux. Similarly, a 't Hooft loop introduces a test magnetic charge along 
a given path, and measures the response to electric flux. Their behavior is dual to one 
another. At zero temperature in a SU(A^) gauge theory without quarks, the condensation of 
magnetic charges confines electric charge, so the Wilson loop has area behavior, and the 't 
Hooft loop, perimeter. Conversely, at temperatures above that for deconfinement, magnetic 
charges do not condense, and electric charge is not confined; hence a (thermal) Wilson loop 
has perimeter behavior, and the (spatial) 't Hooft loop, area [5]. 

The 't Hooft loop does not have a simple representation in terms of the vector potential for 
the gauge field, but in the deconfined phase, the area behavior of the spatial 't Hooft loop can 
be simply understood [5] as a Z(A^) interface [6, 7, 8, 9, 10, 11]. Without dynamical quarks, 
a SU(A^) gauge theory has equivalent vacua, which differ by global Z(A^) transformations 
from one another [3, 4]. To probe this, take a box which is long in one spatial direction, say 
of length L in z, and let the two ends of the box differ by a Z(A^) transformation. Thus at one 
end of the box, z = 0, the Wilson line in the imaginary time direction, V exp{ig J AqcIt) 
equals the unit matrix; at the other end of the box, z = +L, the Wilson line is a Z(A^) 
phase, exp(27rz/A^), times the unit matrix. These boundary conditions can be imposed by 
introducing a background field for the timelike component of the vector potential [6, 7, 8, 
9, 10, 11, 12, 13, 14, 15, 16]: 

Ao{z) = Q{z) , (1) 

where Q is a diagonal matrix in color space, and g is the gauge coupling constant. The matrix 
Q{z) is then chosen to vary so that a Z(A^) interface, centered at z = +L/2, forms. Only the 
ends of the box represent allowable vacua, so a nonzero color electric field is generated along 
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the 2;-direction, Ez ~ dzQ{z), and the configuration has nonzero action. By construction, 
the interface is independent of the x and y directions, and so the action is proportional to 
the transverse area. This Z(A^) interface is then equivalent to a spatial 't Hooft loop, in the 
plane of x and y, at z = +L/2 [5]. 

The action for the Z(A^) interface can be computed in weak coupling, and reduces to a 
tunneling problem in one dimension [6]. There is no potential for Q classically, but one is 
generated at one-loop order, and so the action for the associated instanton is not ~ 1/(7^, 
but ~ \l [6, 17]. A derivative expansion can be used to compute, because along the 
^-direction, the width of the interface is proportional to the inverse Debye mass, ~ 1/((7T), 
which is large relative to the typical thermal correlation length for massless fields, ~ l/(27rT). 
Consequently, at the outset one computes for a field Qiz) which is constant in as effects 
from the variation in z enter through corrections which are of higher order in g. Corrections 
to the interface tension have been carried out to ~ [7], and are underway to ~ g^ [8]. 
This is to be compared with the free energy (where (5 = 0), which has been computed to 
~ / [18, 19, 20, 21, 22]. 

The interface tension is measurable through numerical simulations on the lattice [10] . This 
includes simulations which model the behavior in real time [11]. In this paper we address 
a related problem semiclassically, by computing the simplest real time response functions, 
for a spatial 't Hooft loop, in weak coupling. To do this, we expand to quadratic order 
about the background field in Eq. (1), taking Q to be a constant matrix, and analytically 
continue to real time. As for quantities computed in equilibrium, near equilibrium there is 
a natural division between momenta which are "soft" , with components of order the Debye 
mass, ~ gT ^ and "hard", ~ T. The simplest real time response functions are the quark 
and gluon self-energies, computed in the hard thermal loop approximation, for soft external 
momenta [23, 24, 25]. 

We perform this computation in order to develop techniques which will enable us to 
address a problem of much broader interest. Resummations of perturbation theory appear 
to break down by temperatures several times the critical temperature [20, 21]. The obvious 
guess is to assume that since ~ 150 — 200 MeV in a gauge theory, that the theory 
has entered a nonperturbative regime by this point, where the QCD coupling is large [2]. 
While the former must be true — confinement cannot be seen in perturbation theory — 
computations of an effective theory for the pressure find that the coupling is moderate even 
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at Tc, with a'fiTc) = g'^^iTc) / (in) ~ 0.3 [22]. This is because in imaginary time, the typical 
"energies" are large, multiples of 2'kT [18]. 

Why, then, does deconfinement occur at moderate coupling? Deconfinement is an order- 
ing of global Z(A^) spins, and is measured by the trace of the thermal Wilson line, which is 
the Polyakov loop [3, 4]. In the fundamental representation, without quarks the expectation 
value of the Polyakov loop vanishes below Tc, and approaches one at high temperature. 
Since the Polyakov loop is not equal to one whenever Q 0, one way to model the re- 
gion where the Polyakov loop is not near one is to assume a nontrivial distribution of Q's 
[11, 26, 27, 28, 29, 30, 31], which in Ref. [32] we term a "semi" Quark-Gluon Plasma (semi- 
QGP). At least at infinite A^, it is easy to model the confined phase, as a distribution which 
is fiat in Q. This implies that the expectation value of the Polyakov loop in any nontrivial 
representation vanishes, whether or not the loop carries Z(A^) charge. This is also consistent 
with how the 't Hooft loop must change near T^. At high temperature, Z(A^) interfaces are 
rare, and in infinite volume the theory lies in one Z(A^) domain. As T — * T+, though, the 
interface tension decreases, Z(A^) domains become plentiful, and the Z(A^) spins are disor- 
dered. The decrease of the Z(A^) interface tension near Tc has been confirmed on the lattice 
[10]. 

A semi-QGP can be shown to occur in one unphysical limit. Let the spatial volume be a 
sphere of hadronic dimensions, so small that by asymptotic freedom, the coupling constant 
runs to a very small value [33]. Although systems at finite volume cannot have phase 
transitions, they can if the number of colors is infinite. If R is the radius of the sphere, then 
even when g"^ = 0, aX infinite there is a deconfining phase transition, of first order, when 
Tc = c/R, where c = l/log(2 + -\/3) [33]. At T+, the expectation value of the (renormalized) 
Polyakov loop, in the fundamental representation, is exactly 1/2. The Q-distribution for 
the constant mode on the sphere reduces to a type of matrix model, which can be computed 
analytically about Tc [33]. Further, since all resonances are of zero width at infinite A^, the 
Hagedorn temperature is a precise quantity. On a small sphere at A^ = cxd, at zero coupling 
the Hagedorn temperature coincides with Tc. Perturbative corrections move Tc below the 
Hagedorn temperature, and the loop at T+ away from 1/2, by an amount ~ {g'^{R) N)"^ [34]. 
Presumably, the expectation value of the (renormalized) Polyakov loop, in the fundamental 
representation, goes from 1/2 at T^, to near one, by temperatures which are a few times 
Tc. This is then the semi-QGP on a small sphere. 
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Of course this might be an artifact of working on a small sphere at infinite A^. In a 
large spatial volume, one must look to numerical simulations on a lattice. In any volume, 
the Wilson line, and so the Polyakov loop has ultraviolet divergences, so that the bare loop 
vanishes in the continuum limit. A nonzero value in the continuum limit is obtained after a 
type of mass renormalization [27, 35, 36, 37, 38, 39, 40]. On a small sphere, when g^{R) <^ 1, 
one can renormalize the loop perturbatively. In a large volume, on the lattice there are two 
methods of renormalizing the Polyakov loop [27, 35], which now agree up to the numerical 
accuracy [40]. The most precise measurements are for a SU(3) gauge theory without quarks 
[40]. From Fig. 1 of [40], in the triplet representation the expectation of the renormalized 
Polyakov loop vanishes below Tc, and is ~ 0.5 at Tc. It then rises to it rises rapidly, and 
is ~ 0.9 by 2.0 Tc. It then rises slowly, reaching ~ 1.1 by ~ 4Tc. From 4 to 12Tc, its 
expectation value is flat. This suggests that there is a nontrivial Q-distribution about Tc, 
which is relevant up to temperatures which are 2 — 4Tc. Above 4Tc, any nonperturbative 
effects from the Q-distribution appear negligible, consistent with the success of resummations 
of perturbation theory [21, 22]. In a SU(3) gauge theory with quarks, the simulations are of 
more limited accuracy, but a similar picture emerges [39]. The principal difference is that 
the expectation value of the renormalized triplet loop is nonzero even below Tc [39]. 

At present, it is not known what the Q-distribution is in a gauge theory, even without 
quarks. This would correlate the pressure with the expectation value of the (renormalized) 
Polyakov loop(s). Such a distribution might be obtained from numerical simulations, in 
both the original and an effective theory [30]. 

Even without knowing the full Q-distribution, though, one can take the flrst steps toward 
the response functions in real time, which we do here. For either a Q-distribution, or a 't 
Hooft loop, one begins by computing the quadratic fluctuations about the background fleld 
in Eq. (1). We then analytically continue the fluctuations, computed in imaginary time, to 
real time [25, 41]. (For a numerical approach to quantities in real time, see [11].) 

We flnd that the quark self-energy is a trivial extension of that for Q = 0: it is equal to 
a thermal mass squared, ~ fi'^T^, times the same hard thermal loop as in zero fleld. The 
gluon self-energy is different, though. Besides the usual hard thermal loop, proportional to 
a thermal mass squared, ~ g^T"^-, there is a new piece, ~ g^T^- The function is similar to 
other hard thermal loops, but is novel. We suggest that it arises because of the background 
color electric fleld in a 't Hooft loop. 
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The paper is organized as follows: In Sec. II, we discuss how the double line notation, 
which is standard at large N [42] , can also be used easily at finite N [43] . Our discussion is 
elementary, but is absolutely essential to being able to compute in an arbitrary background 
field with Ao ^ 0, Eq. (1). 

Section III introduces the background field calculation for constant Aq. The perturbative 
rules in the 't Hooft basis are given in Sec. Ill A. These are nothing more involved than 
the usual perturbative rules, with a simple "shift" in the energies, ~ Q. Sec. IIIB gives 
expressions, useful for calculation, in terms of a "mixed" representation, working with spatial 
momenta and Euclidean time. In Sec. Ill C, we follow Furuuchi [41] and discuss how to obtain 
scattering amplitudes, in real time, from those computed in imaginary time. 

The computation of the hard thermal loop in the quark self-energy is given in Sec. IV. 
We go through this example in some detail to develop familiarity with computing diagrams 
when Q ^ 0. The hard thermal loop in the gluon self-energy is computed in Sec. V. In 
Sec. V A we discuss the Q-dependence of the hard thermal loops in two tadpole diagrams, 
which are independent of the external momenta. In Sec. V B we consider the hard thermal 
loops which arise from diagrams which exhibit with Landau damping. With these examples 
in hand, computing the one point gluon function, in Sec. V C, the quark contribution to the 
gluon self-energy. Sec. VD, and that of ghosts and gluons to the gluon self-energy. Sec. VE, 
is relatively straightforward. 

In an appendix we draw some distinctions on the differences between Z(A^) and U(l) 
interfaces [29]. 

II. THE DOUBLE LINE NOTATION AT FINITE iV. 

In order to compute efficiently, it is useful to have a convenient basis for the generators of 
SU(A^). In this section we follow Cvitanovic [43] to show how the usual double line notation, 
which is familiar at large N [42], is also natural at small A^. 

Our purpose here is to establish the notation that we need to compute in the presence 
of a constant, background field for Aq, Eq. (1). We note that at A^ = oo, Aharony et al. 
computed the free energy with 7^ to three-loop order, ~ g"^, on a small sphere [33]. At 
finite A^, recently Korthals-Altes used the double line notation to compute the free energy 
for 7^ to ~ (jf^ in supersymmetric theories [8]. 
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The standard choice for the generators of a gauge group, A"^, is to take a complete and 
orthonormal basis, 

tr (A^ A^) = ^ (2) 

A and B refer to adjoint indices, which for SU(A^) run from A, B = 1,2, . . . , {N"^ — 1). 

We denote indices in the fundamental representation by a, 6 . . . = 1, 2, . . . , A^. Indices 
in the adjoint representation are then denoted by a pair of fundamental indices, (ab). The 
basic quantity which we need is a projection operator which ties together upper and lower 
adjoint indices: 

n1 = 5^ - ^ 4d . (3) 
For an arbitrary matrix M, adjoint indices are raised and lowered by flipping the order in 
the pair: M""^ = Mf,a. The indices flip because off-diagonal generators are ladder operators. 
With this convention, the projection operators between a pair of upper indices, or a pair of 
lower indices, are 

^ab ^ ^ab4c ^ p^^^^ (4) 

It is obvious that Eq. (3) represents a projection operator 

-pff vii = . (5) 

While we have to be careful in lowering and raising adjoint indices, we can raise or lower 
single indices without concern, 5°"^ = — ^ab- The second term in Eq. (3) ensures that it is 
traceless in either of the two pairs of adjoint indices, 

5ab Vf, = 5^" Vf, = . (6) 

While we call (ab) an adjoint index, this terminology is somewhat misleading. SU(A^) 
has A^^ — 1 independent generators, but there are obviously A^^ values for the index (ab). 
By using projection operators (or in the terminology of [43], invariant tensors), the resulting 
basis is overcomplete, with one extra generator. 

While we give explicit expressions for all quantities, it is handy to use a diagramatric 
notation [42, 43]. For SU(A^), lines always carry an arrow, with fields in the fundamental 
representation represented by a single line, and those in the adjoint, by a double line. For an 
adjoint index (ab), we adopt the notation that for upper indices, a is outgoing and b ingoing; 
for lower indices, a is ingoing and b outgoing. This reversal is necessary so that upper and 
lower indices are contracted accordingly. 
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FIG. 1: Projection operator, P"''''^'^ 
a b 



a b a b 

1 ^ 

d c » (i 



FIG. 2: Generator for SU(iV), times \/2. 

Generally, any confusion with indices is dispelled by drawing the corresponding diagram. 
For example, the projection operator of Eq. (3) is illustrated in Fig. 1, and is drawn like a 
gluon propagator. 

With these conventions, the generators of the fundamental representation are just pro- 
jection operators, 

(t^'U = ^ K ■ (7) 

Note that the upper pair, (ab), refers to the index for the adjoint representation, while the 
lower pair, (cd), refers to the components of this matrix in the fundamental representation. 
This is illustrated in Fig. 2. While this is also a projection operator, we draw it differently 
from Fig. 1, distinguishing between the adjoint indices, on top of the diagram, and the 
matrix indices for the fundamental representation, on the two sides. 

As each generator is a projection operator, the trace of two generators is again a projection 
operator: 

tr (t"^ t"^) = - Vff Vt = - V'^'"^ . (8) 
2 2 

We now make an extended comment about the normalization of generators which is implied 
by Eq. (8). While mathematically elementary, at least we found it confusing at first. 

The off-diagonal generators are the customary ladder operators of the Cartan basis. That 
is, for a b, they are normalized as in Eq. (2), i.e., 

tr {t^' &) = ^ ; (9) 

here a and b are fixed indices, with no summation convention. In SU(A^), there are N{N — 1) 
off-diagonal generators. 
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The only difference lies in the choice of the diagonal generators (which is the Cartan 
subalgebra, the space of mutually commuting generators). The Cartan basis includes one 
generator proportional to 

^^''=^ I ^ " I > (10) 




NV2 

where Iat-i is the unit matrix in — 1 dimensions. In the Cartan basis, the corresponding 
matrix is = — a/A^/ (A^ — 1) t^^ , where the overall constant is required so that A^ obeys 
Eq. (2). In the Cartan basis, the other diagonal generators are like t^^, but for smaller A^. 
For example. 



^ l/V-2 ^ 



A^-1 - ^ 



^2{N-l){N-2) 



-(A^-2) 
y Oy 



(11) 



and so on. We denote the A^— 1 diagonal generators in the Cartan basis as A'^, a = 2, 3, . . . , A^. 
While orthonormal, this basis clearly treats the different diagonal elements on an unequal 
footing, with the A^*'^ element occuping a privileged position. 

In contrast, for the double line basis the diagonal generators are just permutations of one 
another: start with Eq. (10), and simply shuffle where the factor of — (A^ — 1) lies along the 
diagonal, e.g.. 



-1 /-(A^-1) 



and so on. The t^^, . . . ,t^^ are a set of A^ diagonal generators, which manifestly do not 
treat any diagonal element different from any other. This is only possible only because they 
are not independent, that is, their sum vanishes, 

N 

^ t"" = . (13) 

a=l 

Consider the example of two colors, where the double line basis has four generators. 
There are two off-diagonal, ladder generators. 





t'=^ ■ (14) 
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FIG. 3: The product of two generators, times two. 
For the diagonal generators, from Eqs. (3) and (7) there are two contributions. 

As expected, t^^ is proportional to the Pauli matrix a^, but the constant appears wrong. 
While the ladder generators are normalized as in Eq. (9), tr(t^^t^^) = 1/2, the diagonal 
generators satisfy tr (t^^)'^ = tr (t^^)^ = 1/4, instead of 1/2; further, tr (t^^t^^) = —1/4. 

This normalization is correct, and arises because the basis is overcomplete. From Eq. (8), 
the trace of a product of two generators is itself a generator. Thus the trace of a given 
diagonal generator, squared, is 

tr in' = ^1 - ^) ' (16) 

with no summation over a; for N = 2, this = 1/4. Further, the trace between two different 
diagonal generators is nonzero: 

tr t"^ = ^) ; (17) 

with no summation over a or b, and a ^ h\ for N = 2, this = —1/4. Thus in this basis, 
the (peculiar) normalization of the diagonal generators arises because the generators are 
projection operators. Of course the diagonal and off-diagonal generators are orthogonal to 
one another in the usual manner. 

In the absence of a background field, that the Cartan basis chooses a preferred direction 
amongst the diagonal generators is of no concern. The Cartan basis is also convenient when 
computing the properties of a Z(A^) interface [6], since then the background field is along 
t'^^ , and treating the N^^ diagonal element as special is natural. 
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In the presence of an arbitrary background field, Eq. (1), though, where Q is a diagonal 
matrix, {Q)ab = Q^'^ab, Eq- (32), the double line basis is more useful. In particular, all 
covariant derivatives are simple. In the fundamental representation, Q acts linearly upon 
fields ip. Then ip is like a column vector, so if ipa is the a^^ element, 

QiJa = Q''iJa. (18) 

In the adjoint representation, the covariant derivative involves a commutator. The commu- 
tator of Q with any generator, though, is just that generator times a difference of Q's: 

[Q, t"''] = {Q" - Q^) = Q^^ . (19) 

This is clear: if a = 6, t°-^ is a diagonal matrix, and the commutator of 1°"°" with another 
diagonal matrix, Q, vanishes. If a 7^ 6, only the first term in (t"'^)cd, ~ ^c^d; Eqs. (3) and 
(7), contributes, to give Eq. (19). We introduce the notation = — Q^, which we shall 
use extensively. 

As a consequence of Eqs. (18) and (19), we find in Sec. Ill A that with the double line 
basis, perturbation theory for Q 7^ is a trivial generalization of that for Q = 0: just a 
constant, albeit color dependent, shift in the energies. Energies which carry color indices 
were introduced when the determinant in a background constant Aq field was first computed, 
in Appendix D of Ref. [12]. 

The double line basis is useful in other ways. While admittedly perverse for two colors, 
for three or more colors it is a very efficient means of deriving various identities amongst 
generators of the gauge group. For example, the product of two generators is 

{t^' t^') ^^ = \ Pff = \ [5: 51 - ^ {5: 5) + 51 53 + ^ S'^' 5^' 5e/) , (20) 

as illustrated in Fig. 3. This certainly shows how writing down all of the indices is more 
tedious than just drawing the corresponding diagram. By tying the sides of the diagram 
together, representing summation over the matrix indices, we obtain the normalization con- 
ditions above, Eqs. (16) and (17). By tying the adjoint indices together on top of the 
diagram, we also obtain the familiar identity: 

E r '"L = ^ ■ (21) 

a,b=l 
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FIG. 4: Structure constant for SU(A^), times —i\f2. 

In principle, we really should contract the adjoint indices with a projection operator. Since 
the generators are traceless, though, Eq. (13), the projection operator reduces to ordinary 
Kronecker deltas. 

The product of three arbitrary generators can be written out, similar to Fig. 3. By 
drawing diagrams, it is easy obtaining the standard relation, 

^ . . . 1 . 

(22) 



^ 1 
2N 



c,d=l 

where we leave the matrix indices implicit. 

More useful is to take a trace of the product of three generators. The antisymmetric 
combination is proportional to the structure constant of the group, 

N 

(23) 



\j^ab ^cdj _ ^ j:{ab,cd,ef) ^fe 

ej=l 



and is simple. 



j!(ab,cd,ef) _ _J_ ^^ad ^cf ^eb _ ^af ^cb ^ed^ 
a/2 



as illustrated in Fig. 4. In this basis the structure constants satisfy the identity: 



N 



j{a.b,ef,gh) j(cdje,hg) _ jy -p' 

e,/,S,/i=l 



ab,cd 



(24) 



(25) 



It is also easy to draw the diagram for the symmetric structure constant, 

^{abM,ef) _2 tr {e^f^.e^}) 



as illustrated in Fig. 5. 
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FIG. 5: Symmetric structure constant, times -\/2- 



For higher representations of the group, instead of diagrams with just two hues, one ob- 
tains diagrams with many hues, or "birdtracks" . For a careful discussion of the classification 
of arbitrary representations of Lie groups by means of birdtrack diagrams, see Cvitanovic 



III. COMPUTING IN BACKGROUND FIELD GAUGE 
A. Propagators in a background field 

In this section we develop the perturbative rules in the appropriate background field [16]. 
At tree level the Lagrangian is 



We assume there are Nj flavors of quarks, ip, in the fundamental representation of the gauge 
group; the covariant derivative in that representation is = — igA^. The field strength 
tensor G^^ = [D^, D^]/{—ig) = d^A^ — d^A^ — ig[Af^,A^]. The covariant derivative in the 
adjoint representation is = — ig[Afj_, . ]. We work in Euclidean spacetime, with a 
positive metric. The gamma matrix is Hermitian and satisfies {'Jfi,'yu} = 26^^. 
We expand about a background field A^', 



[43]. 



(27) 



A, = A';l + B, 



(28) 
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where denotes the fluctuation. The classical covariant derivative is then D'^^ = d^ — igA'^^, 
etc. The gauge fixing and ghost terms are chosen to be those for background field gauge, 
with gauge fixing parameter ^: 



C„ = ^ tr {Dp,Y - 2 tr {vD^^D.v) , (29) 



with rj the ghost field. 



The inverse propagators follow directly. That for the quark field is p'^^ + m. The inverse 
propagator for the ghost is — (Z^^')^, while that for gluon fluctuations, 5^, is 

= - {DfY V + - ] D'^iDt + ig [G% . ] 

Most of our calculations are done assuming a background field which is constant in space- 
time. Notice, however, that the last term in the inverse gluon propagator is proportional to 
the field strength tensor of the background field. This will be important in understanding 
novel terms for gluon hard thermal loops in the presence of an interface. 

The covariant derivative in the fundamental representation enters into the quark inverse 
propagator, while that in the adjoint representation enters into the ghost and gluon inverse 
propagators. We now compute at a nonzero temperature T in the imaginary time formalism, 
where the Euclidean time r : — 1/T. The energies are then 

Po = 2n7rT, bosons ; po = {2n + l)7rT, fermions . (31) 

We use a tilde for the energies and momenta of fermions, to distinguish them from bosons. 

We take the background field as a constant, diagonal matrix for the timelike component 
of the vector potential: 

At = -^Q ; (Q)ab = Q"5,,; (32) 
as an SU(A^) matrix, the sum of the Q^'s vanishes, 

N 



^g'^ = 0. (33) 



a=l 



As discussed in Sec. 11, the great virtue of the double line notation is that the covariant 
derivatives in a field like Eq. (32) are trivial. For fields in the fundamental representation. 
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pa _ S""^ 

b * CL -if"' + m' 

pab 



c < b 

P~ ab 

^ * 1"^"^, r"" = f V - (1 - 1^ I TT^L^ v^'^^" 

c * b 



Id-^^a 

NrJ^h 



1^'^ ''^ (^pab^2 j (^paby 

FIG. 6: Quark, ghost, and gluon propagators. 



the covariant derivative in the background field is D'^^ipa = ~iP^'4^a , Eq. (18). If this 
covariant derivative acts upon a quark field, is a momenta with one color index, 

^; = (po + Q^p). (34) 

Since the background field shifts the Euclidean energies, it is convenient to write 

P; = P^ + Q\ = Q'' + 7rT . (35) 

That is, we treat all momenta as bosonic, which we can easily do by just putting the change 
in the boundary condition for fermions, versus bosons, into part of the background field. 
The covariant derivative acts upon fields in the adjoint representation as D'^^t"'^ = 

_^pab^ab ^ ^igy ^OT boSOUlc ficlds, 

P^' = {po + Q^-Q\p) = {Pl\ P) , (36) 

and involves an adjoint color index, {ah). 
To be explicit, the quark propagator is 

( r{P) ^\-P) ) = —~ ; (37) 



the ghost propagator. 



and the gluon propagator. 



( r^»^(P) rJ^'i-P) ) = ^ V^''^' ; (38) 



pab pab \ i 

( BfiP) Bfi-P) )= I 5,. - (1 - -j^j ^ V^'^'' , (39) 
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These are illustrated in Fig. 6. 

There are several matters of notation to attend to. All of these sound more complicated 
than is true after drawing the corresponding double line diagram. In Eqs. (34) to (39), we 
adopt the convention that color indices shared between momenta and projection operators 
are not summed over. 

Implicitly, a quark line carries two arrows: one as a Dirac particle, and one for color. 
Either Pq = Po + Q"', if the directions coincide, or (P")o = Po — Q"^, if they do not. As 
long as one is consistent with directions, though, this does not really matter. For the quark 
propagator, this leads to the obvious identity: 

= . (40) 

-i^'^ + m -iJf^ + m 

Similarly, for fields in the adjoint representation, where P^^ = Po + Q°' — Q^-, we define the 
left index so that it is in the direction of the momentum. Thus if we change the direction, 
P' = —P, then we must also reverse the order of indices, and (P')*"^ = —P"-^. For the ghost 
propagator, for example, 

^ 'pabjCd ^ 'pab,cd (41) 



The same is true for the gluon propagator. 

The vertices between quantum fluctuations, and if), are also simple, taking care of 
indices and the like. The vertex between a quark, antiquark, and a gluon is obtained by 
taking the derivative of the action, S = J d'^xC, as 

— =a = ^9 it"') , la . (42) 

The order of the gluon indices flip, from the left to the right hand side, because B^ = t'^'^B'^f. 
The vertex between a ghost, antighost, and a gluon is 

ig (p-b)^ . (43) 



6r]f^{R) 6Bff{Q) 5ff''{P) 
The three gluon vertex is 

_ 5S _ 

~ 5B{\R) 5Bi'^{Q) 5B^-{P) ~ 

where 



g f^''''"'''^^ P^.uxiP''', R'^) , (44) 



r,.x{P''\ Q"'. R'^) = {P^' - ) 6,. + (Q^' - R"/) 5.x + (K^ - ) ■ (45) 
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These are all the usual vertices, with the replacement of ordinary momenta by momenta 
which carry color indices. Again, in Eqs. (43) and (45), color indices shared by momenta 
are not summed over. 

Momenta with colored indices satisfy momentum conservation as usual, so in Eq. (45), 

+ Qf + < = . (46) 

Now consider one of the external momenta, say , contracted with the three gluon vertex. 
This satisfies the identity: 

Rl^ r,,,{P^\ Q^", Wf) = A-liQ^') - A-l{P'^') , (47) 

where A^^ is the transverse piece of the inverse gluon propagator: 

A-l{P)=5,^P'-P,P,. (48) 

This is the same identity as for Q = 0, Eq. (2.50) of Ref. [24]. This is useful because as for 
(5 = [24], it can be used to show that hard thermal loops are independent of the gauge 
fixing condition. 

The four gluon vertex has the usual form, a sum over products of structure constants: 

1^1 = _ ^2 / f{ab,cd,ij) r{ef,gh,ji) / e r _ r r N 

6B'.^{S) 5Bi\R) 5Bi^{Q) 5B^;{P) ^ ^ ^ '"""'^ 

+ f(-b,9h,ij)ficd,em ^^^^^^^^ _ _ 

(49) 

We conclude by noting that Eq. (41) can be used to simplify insertions of ghost or gluon 
lines in loop diagrams. If a gluon ties onto a quark line, then in Feynman gauge this enters 
as 

'r>ab,cd 

(p^2 ' (50) 
where we neglect the rest of the diagram. The projection operator in the gluon propagator 
is a sum of two terms, Eq. (3). Since the gluon appears inside the loop, the c and d indices 
are summed over. Unlike the a and b indices, which also enter through P"'', this is the only 
place where c and d indices enter. Since the generators are traceless, though, Eq. (13), any 
contribution from the second term in the projection operator, —5"'' vanishes. Hence 
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in the gluon line, we can replace the projection operator by the first term, which is just a 
Kronecker delta, 5"''^5^^, so Eq. (50) becomes 



The same is true in any gauge. It is also true for a gluon tied to either a three gluon 
or four gluon vertex, since in each case, Eqs. (45) and (49), the projection operator in 
the gluon propagator ties onto a factor ^ and that any such structure constant 

involves a commutator of t^'^, Eq. (23). The same holds for a ghost propagator tied onto a 
ghost antighost gluon vertex. This useful simplification was first seen in Eq. (21) in Sec. II: 
generators can be contracted not with projection operators, but just with ordinary Kronecker 
deltas. 

B. Propagators in a mixed representation 

In this section we discuss a useful trick for computing scattering amplitudes, starting in 
the imaginary time. To avoid unnecessary complication, we replace the color matrices, either 
Q"^ or Q"''', by a single background field, Q. This is identical to considering the propagation 
of an electron in QCD, in the presence of a background field Aq ~ Q/e. The extension to 
QCD is automatic, as will be clear from the examples which follow in later sections. 

We then introduce a "mixed" representation for the propagators [23, 24, 25]. For the 
spatial directions, one works as usual in momentum space, but for the time direction, instead 
one stays in coordinate space. For example, consider a propagator. 



energy po = ^imT, for integral n. For the time being, we also assume that < r < 1/T. 
The sum is performed by contour integration in the complex po plane. There are two poles, 
for = —Q±iE, which give 




(52) 



E is the energy, typically E = a/p" 



2 + m^. We assume the field is bosonic, so the Euclidean 




(53) 



Here n{E) is the usual Bose- Einstein statistical distribution function. 



(54) 
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so the only change for Q 7^ is the change in the statistical distribution function, 

Notice that this is the only place where Q enters into Eq. (53): the propagators in Euclidean 
time, ~ exp(=i=£'r)/2i?, are identical to that for Q = 0. 

This can be understood by recognizing that the parameter iQ enters exactly like a chem- 
ical potential, albeit one which is imaginary. Because of this, the associated statistical 
distribution functions, n{E =]= iQ)^ are complex valued. It helps to rewrite Eq. (53) in a less 
compact form, 

Z\Q(r, E) = ^ ((1 + n{E - zQ)) e"^- + n{E + iQ) e+^^) . (56) 

This form is physically more transparent. The first term, with propagator e"^^/2i?, is 
proportional to 1 + n{E — iQ). The 1 is the contribution in vacuum, while n{E — iQ) 
represents the induced emission of a particle, with energy E and chemical potential +iQ, 
into the thermal bath. The second term, with propagator exp(+i?r)/2i?, is proportional to 
~ n{E + iQ). This represents absorption of a field with energy E., and chemical potential 
for the antiparticle, —iQ, from the thermal bath. 

This expression for the propagator is the same as for a real chemical potential, /i, except 
that /i is replaced by iQ. As when /i 7^ 0, the statistical distribution functions are modified, 
but the energies of the system remain unchanged. The same is then true for an imagi- 
nary chemical potential, iQ. This is why the form of the propagators in imaginary time, 
exp(±i?r)/2£', are unaffected by Q. We argue in the next subsection that this remains valid 
for propagation in real time as well. 

We make some remarks to help illuminate the meaning of the statistical distribution 
functions when < Q < vrT. First, when (5 = 0, the Bose-Einstein distribution function is 
singular as — 0, n{E) ~ T/E. This singularity is related to the phenomenon of Bose- 
Einstein condensation at low temperature. In contrast, whenever Q 7^ 0, the distribution 
function is regular as the energy vanishes, n{{] — iQ) = l/{exp{—iQ/T) — 1). This includes 
ordinary fermions, when Q = vrT, and n(— ivrT) = —1/2. 

It is also helpful to consider adding a real chemical potential, fi, in addition to iQ. For 
ordinary fermions, Q = vrT, the Fermi-Dirac distribution function with /x 7^ is 

n{E - /i) = -n{E - inT - /i) = ^^e_^It ^ i ■ 
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In the limit of zero temperature, if £" > /i, n = 0, while if £" < /i, n = 1. For antiparticles, 
n{E + fi) = for any E as T ^ 0. This is just a Fermi sea, and represents a net excess of 
particles over antiparticles. 

When Q ^ ttT, if there is also a real chemical potential, /x, the associated statistical 
distribution function is 

n{E-fi- zQ) = ^^E-^.-l)/T _ 1 • (58) 
Taking Q = 27rTg, as T ^ 0, if > /i, n = 0. However, if E < fi, n = —1; the negative 
sign of n is natural, see Eq. (57). For antiparticles, n{E + fi + iQ), one finds n = for all 
energies. Thus a real chemical potential introduces an asymmetry between particles and 
antiparticles for all Q 7^ 0. 

Speaking loosely, when < Q < ttT particles behave with something like fractional 
statistics. This analogy is not precise, though, merely suggestive. In particular, for both 
cases in which Q arises, a given Q is not physical. For a Z(A^) interface, Q is a function of 
z, and one integrates over all Q{z). In the semi-QGP, there is a distribution of Q's, and it 
is only integrals over the distribution which are physically meaningful. In both cases, after 
summing over all Q's, the usual relationship between spin and statistics is recovered. 

We conclude with some useful identities. The first is 

l+n{E- iQ) = e(-^-^«)/^ n{E - iQ) . (59) 

This is well known for Q = 0, and by construction, must then be true for Q 7^ 0, by simply 
replacing E ^ E — iQ. It applies for either sign of Q. 

The propagator in Eq. (52) is defined for positive t, < r < 1/T. The extension to 
negative values, —1/T < r < 1/T, is 

Aq{t, E) = Y^^{1 + n{sE - iQ sign(r) ) ) e^^l^l . (60) 

s=± 

From this, or directly from Eq. (52), 

Aq{-t,E) = A_q{t,E). (61) 

From this it also follows that 

AQ{r-l/T,E) = e'^l^ AQ{r,E) . (62) 
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This is the generahzation of the Kubo-Martin-Schwinger condition [25] to a background field 
Q. 

In practice we will start with diagrams in momentum space, and then transform a sum 
over pq to an integrals over r's. This is done by using the Fourier transform of Eq. (52), 
which is 

In computation we also require 
Po + Q [^'^ , f d 



(64) 



(pO + Q)2 + ^2 

Integrating by parts, this equals 

_ , (e*(fo+Q)/T^Q(i/r^ ^) _ ^^(0, E)) + y rfre^(f«+«)^ '5:^ ' 

The first term vanishes, since po is bosonic, so e*^°/"^ = 1, and by the condition of Eq. (62), 
for T = 1/T. The second term is easy to evaluate, and gives 



/ dr ((1 + n{E - iQ)) {-tE)e-^^ + n{E + tQ){+iE)e+^^) , 

Jo 

(66) 

Integrals with higher powers of po + Q cire not required, since they can be reduced as 

(Po + Q? , 



{po + Qy + E^ (p, + Qy + E^ 

which can be handled by previous results. 



(67) 



C. Amplitudes in real time 

In this section we follow Furuuchi [41] and discuss how to proceed from amplitudes, 
computed in imaginary time, to scattering amplitudes. 

We remark that Smilga [15] has argued that Z(A^) domain walls are entirely a construc- 
tion valid only in imaginary time, and have no relevance for scattering amplitudes. If Z(iV) 
domain structure is a natural consequence of the system in thermal equilibrium — i.e., in 
imaginary time — then it is difficult to see how they could not be relevant for small fiuctua- 
tions about thermal equilibrium, which is what amplitudes in real time represent. Notably, 
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FIG. 7: Contour in real and imaginary time. 

lattice simulations of the dynamical evolution of Z(A^) domains have been performed [11], 
and as one might expect, are very similar to the evolution in Potts models, to which they are 
closely analogous. We do acknowledge, however, that the lattice time step is not immediately 
related to a physical time. 

Consider the usual contour in the plane of real and imaginary time. The imaginary time 
variable r runs from to 1/T, and represents a thermal ensemble in thermal equilibrium. 
The contour also runs in real time, t, from to oo, and then back again, representing 
fluctuations about the thermal ensemble. This is illustrated in a standard figure. Fig. 7. 
The exact shape of the contour [25] will not matter for our purposes. 

The essential question is on what parts of the contour the background field Aq is nonzero. 
The answer is to take the background Aq field only for the part of the contour in imaginary 
time, and not for the part of the contour in real time. This is absolutely necessary for the 
integrals to be well defined. Since the real time runs from — oo to +oo, and then back again, 
if there was such a background field, it would affect the behavior at large times. 

This is clearest in considering the background Q field as a chemical potential for color 
charge [14]. A real chemical potential alters the initial statistical distribution of the particles: 
for fermions, for example, it represents a net excess of particles over antiparticles, or a Fermi 
sea. While the canonical momenta are shifted by a chemical potential, the evolution in real 
time is by the usual Hamiltonian of the system. Thus the evolution, in real time t, of some 
operator proceeds by the usual Heisenberg relation, = e"'"*-'^*0(O)e~*'^*. 

While the Q field represents an imaginary chemical potential, its effect is only to alter 
the initial color distribution of particles. The canonical momentum is shifted, but not the 
Hamiltonian. 
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In practice, one computes an amplitude with Q dependent momenta: Po for gluons, 
for quarks. The above imphes that one amphtudes are constructed by taking 

pS'=Po + (g"-g'')^-ic^"\ (68) 

for gluons, and similarly for quarks. Here, u) is an energy in real time, and as such, can take 
arbitrary values. 

The division into hard momenta, of order ~ T, and soft momenta, ~ gT, is preserved by 
this procedure. The usual Euclidean po is a multiple of 27rT for bosons. Thus, the shift by 
the fractional amount, ~ Q, does not change this. In contrast, the Minkowski energy o;"^ is 
a continuous variable, and it is consistent to assume that it is soft. 

There is an important subtlety which wc ignore. Usual scattering amplitudes are invariant 
under arbitrary reparametrizations of the fields. It is far from clear that this is true for 
scattering in a fixed Q field. The point is that it is necessary to compute a physical process. 
For a,Z{N) interface, this would be scattering not at a given point, but integrated over the 
entire spatial extent of the interface. Similarly, for the semi-QGP, scattering in a fixed Q 
field could well exhibit unphysical behavior. The physical quantity there are amplitudes in 
which one integrates over the entire distribution of Q's, representing the thermal equilibrium 
state. 

IV. QUARK SELF-ENERGY 

With the formalism in place, we proceed to computing the self-energy for a quark in 
a background field, Q 0. The computations are relatively straightforward, and do not 
exhibit complications which will arise for the gluon self-energy in the next section. We 
go through this example in some detail, so that the reader can develop familiarity with 
computing in the presence of a background field. It also helps to understand the novelty of 
the new terms in the gluon self-energy. 

At one-loop order, the standard diagram is, in our notation, 

- EiP-u = - 9^ (t-'u p^j, (f'u 1 1^ fp::'i^:^^.^. ■ m 

Here K"^ is the momentum of the quark in the loop, and P" — K'^ the momentum of the 
gluon. Thus K'^ is a fermionic momentum, and — K'^ = P" — K'^ is bosonic. 
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FIG. 8: One-loop diagram for the quark self-energy. 
In Eq. (69) the integral is that appropriate for a bosonic field at nonzero temperature, 



Remember that we can take po and /cq fo be bosonic, by using a background field which is 

The color structure reduces immediately. By Eqs. (50) and (51), the gluon projection 
operator can be replaced by an ordinary Kronecker delta 

(71) 

The color structure is illustrated by the diagram of Fig. 8. This is a sum of the planar 
diagram, minus times a diagram in which all indices are equal. 

We wish to extract the hard thermal loop, ~ T^, from Eq. (71). Instead, to simplify the 
discussion, we consider the integral 

in full generality. 

To perform the sum over tt,, it is useful to use the mixed representation. Sec. IIIB. Using 
Eq. (63), we write 

1 /-l/T gi(fco+Q'=)T 



{K 



c\2 



dr ((1 + n{E, - zQ^))e-^'=^ + n{Ek + zQ^)e+^^^j . (73) 



for the quark like propagator, = v fc^, and 

1 fVT ^i^po-ko+Q^'^y 



+ n(Ep_fc + ^Q"^)e+^--'="') . 
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for the gluon propagator, Ep^k = y {p — k)"^- 

The sum over n, where ko = 2TTnT, can then be performed immediately, and gives a deha 
function in time, 



+00 



Tj^e"""^"'"'^ =SiT-T') . (75) 



-00 



We divide this up into four integrals. 



(Pk 1 



All of the terms are integrals over r. The first is 
Jo 

with the other three of a similar form. The integral is easy to do, 

/ Jipo+iQ'"'+iQ''-Ek-Ep_k)/T _ i\ 

Ti{Q',Qn= ]{l + n{Ek-^Q')){l + n{Ep.k-^Qn) ■ (78) 

V ipo + iQ'' -Ek~ Ep^k J 

In the energy denominator, we rewrite 

+ Q- = Q--Q- + Q- + T,T = Q'' . (79) 

For the other terms, though, it is better not to use this. Since Pq is a bosonic momentum, 
exp{ipo/T) = 1, and we group 

^{ipo+iQ'"'+iQ'-Ek-Ep_,,)/T ^ ^-{Ek-iQ')/T ^-{Ej,_k-iQ'"') /T _ ^gg-j 

This is useful because the identity of Eq. (59) can now be brought to bear, so that 

Xi(Q^g"^ 



n{Ek - iQ')n{Ep^k - iQ"') - (1 + n{Ek - iQ'')){l + n{Ep^k - iQ"')) 



ipo + iQ" -Ek- Ep_k (^3;L) 
= ~^ ^ (l + n{Ek-tQ'')+n{Ep.k-tQn) ■ 

This is the exact same expression as for (5 = 0, with the replacement 

Ek Ek — iQ"^ ; Ep_k Ep_k — iQ"""^ . (82) 

This holds for the energy denominator as well: that for Q = 0, ipQ — Ek — Ep_k, becomes ip^ — 
Ek — Ep_k after the shift of Eq. (82), using Eq. (79). Note, however, that this substitution 
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is not universal, and holds only for the X's: in Eq. (76), the residues for the propagators 
remain l/(2i?fc) and l/{2Ep_k), respectively. 

We went through this derivation in detail, because the chemical potentials are imaginary, 
so that some care is in order. Even so, the manipulations for Q 7^ are very similar to those 
for Q = 0. For example, in going from the first line in Eq. (81) to the second line, that terms 
involving two statistical distribution functions, n{Ek — iQ'^)n{Ep^k — iQ"''^), drop out is just 
the usual cancellation between stimulated emission and absorption in a thermal bath. 

The other integrals can be done similarly, and follow from the result for Q = 0, by shifting 
the energies, as in Eq. (82), using the identity of Eq. (79): 

UQ', Qn = -a J ^j, - iQ') - n{Ep_, + zQ-)) , (83) 

ip^ - Ek + Ep^k V J 

X3(Q^ Qn = ^ Ue, + tQ'^) - n{Ep^k - iQn) , (§4) 

tpg + Ek + Ep_k V y 

The extension to the original integral of Eq. (71) is immediate. The term ~ ■ 7 is the 
same as above. That ~ k'^'-y^ is evaluated using the identity of Eq. (66). In this context, this 
tells us to replace — > ^lE^- For the terms with positive energy, Xi and X2, —lE^ enters; 
for those with negative energy, X3 and X4, +iEk. 

We next extract the hard thermal loop from Eq. (71). These are terms where the external 
momenta are soft, p ~ gT, and are as large as the corresponding term at tree level. In the 
quark self-energy, then, the hard thermal loops are ~ g^T"^ /p. 

The loop momenta in a hard thermal loop are hard, k ^ T. In the statistical distribution 
functions we can then approximate Ep_k ^ E^ = k. The dominant terms arise from the 
energy denominators with Landau damping, I2 and X3. These terms are dominant because 
the energy denominators are a difference of large energies, and thus are small: 

ipo" ± {Ek - Ep.k) ~ ^Po ±Pcose , (86) 

where cos 6 = k ■ p. 

Notice that these terms are only dominant if — and only if — the energy ip^ is small 
after analytic continuation. Sec. IIIC. If only ipQ, and not ipg, were small after analytic 
continuation, then iQ"- would be a hard momentum, and there would be no hard thermal 
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loop in the diagram: everything would be a correction which is suppressed by at least ~ g 
relative to the propagator at tree level. 

Introducing the vector K = {i, k), = 1, the hard thermal loop in the quark self-energy 
becomes 

- SiP'^U ^ irnl)ab SEiP'^) , (87) 

HTL 

where the ~ sign indicates that the hard thermal loops are equal, but not (necessarily) 
terms beyond that order, and 

the angular integral is over all directions of the unit vector k. The function SIJ{P) is identical 
to that for Q = 0. In the quark self-energy, this function is multiplied by a thermal quark 
"mass", which is a function of Q: 

c=l 



where 



o i-oo 

A{Q) = -— dk k {n{k - iQ) + n{k + iQ)) . (90) 

We normalize A{Q) in anticipation of the final result. For ordinary hard thermal loops, the 
integrals over the hard, loop momenta decouple into an angular integral times an integral 
over the statistical distribution functions. Equation (87) shows that this remains true for 
the quark hard thermal loop when Q 7^ 0. The same is also true for the gluon hard thermal 
loop. 

For the quark self-energy, in the end the color structure is no different than the propagator 
at tree level, ~ 5ab- New color structures do arise for the hard thermal loop in the gluon 
self-energy. 

The integrals over statistical distribution functions when Q 7^ are not much more 
difficult than for Q = [6]. We write 

Q = 2'KTq . (91) 
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Then 

/ dkk n{k - iQ) ^ — dk k = ^ — 



oo 

2 ^ 

i=i 



= l-6?(l-g)+iCl2(27rg) , (92) 

where Cln(^^) is the Clausen function. As is typical of similar expressions for a Z(A^) interface, 
this is valid only for < g < 1, but the extension to other values is direct. The imaginary 
term Cl2(27rg) cancels in the sum which enters into A, so that 

A{Q) = 1 - 6g(l - q) . (93) 

The hard thermal loop in the quark self-energy is identical to that for Q = 0, up to the 
change in the thermal quark mass. To check that one obtains the usual value for Q = 0, 
remember that when all g = 0, Q'^ = ttT, or 1/2. Since A{jS) = 1 and A{l/2) = -1/2, 

_ 1 o2m2 

mlS)a, = V • ^^^^ 

There is a simple interpretation of the function A{Q). Although classically there is no 
potential for Q, in the presence of a background field Q, a potential is generated at one-loop 
order. For convenience we normalize this potential as 

V(g) = \ q\l - qf . (95) 

Then 

AQ) = ^ v(g) . (96) 

Hence the thermal quark mass is naturally the second derivative of a potential, as one would 
expect. 

V. GLUON SELF-ENERGY 

A. Gluonic hard thermal loops: tadpoles 

From the example of the quark self-energy, one might expect that the hard thermal loop 
in the gluon self-energy when Q 7^ is just like that when (5 = 0, with the same functional 
form, and the only change a relatively trivial redefinition of the thermal mass. 
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We shall see that this is not true: there are new terms which arise uniquely for non- 
Abelian gauge fields. These are not present either for fermions, coupled to either Abelian 
or non-Abelian gauge fields, nor for Abelian gauge fields. 

Before delving into the details of the computation in a non-Abelian gauge theory, in the 
next two subsections we discuss the differences between hard thermal loops in the gluonic 
self-energy when Q ^ 0, versus Q = 0. Hopefully this will make the origin of the new terms 
less obscure. 

Hard thermal loops are one-loop diagrams which are as large as the corresponding terms 
at tree level when the external momenta are soft, P ~ gT. For the gluon self-energy in zero 
field, at tree level the inverse propagator is ~ ~ {qTY for soft P. Thus the hard thermal 
loops are diagrams which are ~ g'^T'^ times a dimensionless function of the energy, ip^ = u, 
divided by the spatial momentum, p. 

The simplest hard thermal loop is present in (massless) scalar field theories, and is just 
a tadpole integral, 

d^K 1 

(97) 



(27r)4 12 
Consider the extension of this integral to Q ^ 0: 

f 1 

We use the representation of the propagator in Eq. (66). The sum over k^ generates 5{t), so 
the r integral is trivial, and we are left with a single integral over k = so that Eq. (98) 

becomes 

1 /"OO rri2 

— dk k {n{k - iQ) + n{k + iQ)) = — A{Q) . (99) 
4vr^ Jo 12 

By comparison with Eq. (90), this is the function A{Q) we found for the quark self-energy, 
Eq. (93). (An ultraviolet divergent term at zero temperature, f k dk has been dropped. 
For a scalar field, this is part of a mass divergence. For a gauge field, the sum of all such 
terms vanishes at zero temperature by gauge invariance.) 
A less trivial example is given by the integral 

d'^K ko + Q ttT^ 

As we shall see in Sec. VC, there is a contribution to the one point function of the gluon 
from such a term, where k^ + Q arises from the three gluon vertex. We stress that such a 
diagram does not arise for the two point gluon function. 
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(p^^ + Q _ = ^AW). (100) 



This integral vanishes when Q = because it is odd in ko, and one sums over both 
negative and positive values of /cq- Using the same tricks as above, 

Ao{Q) = -—^ dkk\n{k-iQ)-n{k + iQ)) . (101) 

A term ~ j d^k aX zero temperature is dropped, and certainly vanishes when all such terms 
are summed together. Expanding as in Eq. (92), with Q = 27rTg, 

3 °° 1 

MQ) = ^J2-. sm(27rgj) = g(l - g)(l - 2q) . (102) 
We recognize this as the first derivative of the potential in Eq. (95), 

A(Q) = ^V{Q) . (103) 

dq 

Notice that this term manifestly vanishes when Q = q = 0. It also vanishes for g = 1, 
because this is a vacuum equivalent to g = 0, and for g = 1/2, because this is an extremal 
point of the potential. 



B. Gluonic hard thermal loops: Landau damping 

We next turn to the nontrivial hard thermal loops, which have discontinuities, as a 
function of the external momenta. These arise from energy denominators with Landau 
damping. 

We start by deriving, briefly, how such hard thermal loops arise for the gluon (or photon) 
self-energy when Q = 0. Consider, as a prototype, the integral 

r d^K k^y 

J (2^ K^{p-KY ■ ^ ^ 

For the quark contribution to the gluon self-energy, the factors fc* and y arise from the 
quark propagators. For the gluon contribution, these momenta arise from the momentum 
dependence of the three gluon vertices. There are, of course, other contributions, with 
momenta (^o)^ and k^y, but these can be treated similarly. We take the momenta in the 
numerator to be ~ because the loop momenta for hard thermal loops are hard momenta, 
~ T, and dominate over soft loop momenta, ~ gT. 

This integral is done as for Eq. (72) in Sec. IV. There are two terms which persist at zero 
temperature, Ii and X4, with energy denominators ipQ =F {Ek + Ep_k). If /c ~ T is hard, and 

30 



p ~ gT is soft, then these energy denominators are hard, ~ 2k. Such terms are independent 
of po and p, so that Ti + X4 ^ 2n{k)/k, which just produce a tadpole term as in Eq. (97): 

— / — / dk k^ — ^ = —- . 105 

87r2 J Att Jo k 72 ^ ^ 

Instead, concentrate on the terms I2 and I3, which arise from the denominators with 

Landau damping: 

(106) 

For the residues of the propagators we have taken Ep^k ~ k, but this approximation cannot 
be made so cavaherly in the rest of the expression. Indeed, notice that without any cancella- 
tion, these terms are nominally larger than we expect, ~ j dk k"^ n{k) ~ T^, and not ~ T^. 
However, what enters into Eq. (106) is only the difference of the statistical distribution 
functions. For k T and p ~ (/T, 

n{Ep^k) - n{E,) ^ ^ n{k) (1 + n{k)) . (107) 

Because of this cancellation in the statistical distribution functions, the diagram is not ~ T^, 
but only ~ T^. It is a product of an integral over k, 

1 

dk k"^ n{k) {1 + n{k)) , (108) 



and an angular integral. 



^ ^ ip-k) I — - — . ^ — . I . (109) 

ipo/p — p-k ipo/p + p-k 



2 / Att 



This is the usual hard thermal loop. The integral over k generates the thermal mass for the 
gluon. This is multiplied times an angular integral, which generates a dimensionless function 
of ipo/p- This function has discontinuities on the light cone from the Landau damping of 
massless particles. 

We write the angular integral in Eq. (109) as we do to emphasize its behavior as a function 

of 

X = cos9 = p ■ k . (110) 

The integral over x is from —1 to +1, so a nonzero result must be even in x. In Eq. (109), 
this happens because the difference of statistical distribution functions, and the difference 
of energy denominators, are each odd in x, so the product is even. 
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Now consider the analogous integral when Q ^ 0. As a typical example, consider 

J (p, Q^^Q^)-ij iK+Q,np-K+Q,r ■ 

Here ko and po are both taken to be bosonic momenta, while Qi and Q2 are arbitrary 
background fields. The numerator is chosen to (2A;* — p*) {2^ — p^), which differs terms 
proportional to and p^ from Eq. (104). We need this term to keep a symmetry, Qi ^ Q2. 
This symmetry can be checked by shifting of integration variables, K ^ P — K, since both 
ko and po are, by assumption, bosonic momenta. In the numerator we take ~ k'^y as a 
term which generates the largest terms for hard k; terms ~ k^p^ are down by ~ (7 to this 
term, those ~ p''p' by ^ g^. We keep the numerator of order p because the existence of the 
background field makes the new leading order, which we will see in the following, of order 
g'^T^/p. If one neglects the terms proportional to and p' in the numerator of Eq. (Ill), 
one find a new term which is odd under Qi ^ Q2- To see this symmetry in another way, 
note that the relevant external momentum is P^^ = P + Qi + Q2. Then the hard thermal 
loops in Eq. (Ill) can be rewritten as 

1 r d^K {2k' -p') {2y -p)) _ 1 r d^K {2k' ~p') {2^ ~ f) 

4 J (27r)4 (i^l)2(pl2 _ ;^l)2 - 4 y (27r)4 (i^2)2(pl2 _ ;^2)2 ' (^^2) 

where = K + Qi and = K + Q2. This is not a Q-dependent shift of momenta, which 
would be invalid for arbitrary Q. Instead, P^^ — is just a shift of purely bosonic 

momenta, K + Qi ^ P — K -\- Qi. 

In Eq. (Ill) there are terms which persist at zero temperature, Ii and X4, with energy 
denominators ipl^ =F {Ek + Ep^k)- These are really no different than for Q = 0, Eq. (105); 
the dependence on the external momenta drops out, and one is left with integrals 

/^3 ^ f^i y ]^ 
(2^ Jk^ 2k ^""^^ ~ '^'^ ^ ""^^ ^ '^'^ ^ ""^^ ~ '^'^ ^ ""^^ ^ '^'''^ 

T^ij ^^^^^ 
{A{Q,) + A{Q2)) , 



144 

where we have dropped on the numerator, which is higher order of hard thermal loop 
approximation. This is just like the integral of Eq. (98). The angular integral is trivial, 
~ 6'^ , and the Q-dependence is only through terms ~ T'^A{Q). (As noted before, tadpole 
integrals such as Eq. (100) arise only for the one point gluon function, and not in the gluon 
propagator. This is clear just on dimensional grounds, as that diagram has dimensions of 
(mass)^.) 
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We then turn to the terms with energy denominators which correspond to Landau damp- 
ing, T2 and T3. This part of the integral is 

The statistical distribution functions, and energy denominators, which represent Landau 
damping are the generalization of that for Q = 0, Eq. (106), to Q 7^ 0. These are just 
modifications of Eqs. (83) and (84), replacing the Q's there by Qi and Q2: 



^ ^ _n{Ek-iQi) -n{Ep_k + iQ2) , n{Ep_k - iQ2) - n{Ek + iQi) 
^ + ^ ~ ^ p ^ F + ,-„i2 , w w ' ^^^^> 



Po^ = + Qi + Q2- Now we symmetrize each term with respect to then interchange of Qi 
and Q2- 



1 n{Ek - iQi) - n{Ep_ k + iQi) + ri{Ek - iQ2) ~ n{Ep_k + ^^2) 

ipl^ ~Ek + Ep. 



^2-0 7ri2 — 7^ — ' l-'--'-"^ 



_ 1 n{Ep_k - iQi) - n{Ek + iQi) + n{Ep^k - iQi) - n{Ek + 2(^2) 
' " 2 ipf + Eu-Ep.u ■ ^ ^ 

After symmetrization it is then easy to pick out both the leading, and next to leading 

terms from (116) and (117). 

The leading term is easy. In every term, we approximate -Ep-fc ^ E^ = k, and neglect 
and p' in the numerator, so we find that J'^^{P,Qi,Q2) factorizes into a product of an 

integral over J dk, and an angular integral. The former is 
1 
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dk k'^ {n{k — iQi) — n{k + iQi) + n{k — iQ2) — n{k + iQ2)) 

""^^ (A(Qi) + A(Q2)) , (118) 



12 

using the function ^0 of Eq. (102). In all 



T-^{P.Q,M = '^ {Ao{Qr)+Ao{Q2)) I ^k^kn . / . . + . , ^ . , 
12p J 4:TC yiplyp-p-k ipl^/p + p-k 

(119) 

In the angular integral, produce terms ~ 1 and x^, where x = p ■ k, Eq. (110). The 
angular integral is manifestly even in x, and so does not vanish. This expression can be 
rewritten in a form similar to that of Eq. (87), 

J''{P,Qi,Q2) ^ -r iAo{Qi) + Ao{Q2))sr^{P'') + ... . (120) 

D 
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We introduce the function 



. , f I ; (121) 



An 



\p-k 



K = {i,k), with dQ the integral over k. The hard thermal loop is times J'^^ . Since 
6r ~ 1/P, this is ~ g'^T^/p times a dimensionless function of ipo^/p- This is rather different 
from the hard thermal loops for Q = 0, which are ~ g^T"^ times a dimensionless function of 

ipo/p. 

The origin for this difference is natural, when one considers the propagator in a back- 
ground field, Eq. (30). There are the ordinary terms, ~ (D^^)^, which becomes {P'^^Y in 
momentum space. In addition, however, there is also a term ~ 2ig[G'f^j^, ]. For a Z(A^) 



gGtr^gdA^^]-Td.q{z)r^gT\ (122) 



interface, 

g Gt ^ gdA- 

9 

where we use the fact that the typical spatial momenta for a Z(A^) interface is small, ~ gT. 
Thus this term is larger by ~ 1/(7 than the ordinary terms, g^T"^- 

The hard thermal loop in Eq. (119) is ~ g'^T^/p, which for soft p is ~ gT"^. Thus the 
new hard thermal loop can be viewed as a modification of the term in the background field 
propagator for a gluon. Like any other hard thermal loop, it is as large as the term at tree 
level for soft external momenta. It is just that in a background field, this term is larger than 
expected. 

The term ~ 2ig[G'^^j^, .] is special to a non-Abelian gauge field, since it involves the commu- 
tator in group space. It also has no analogy for a fermion field, either Abelian or non-Abelian. 
This explains why it did not appear in previous examples. 

Having such a term is special to computing for a background field with a nonzero color 
field; for a Z(A^) interface, it is a nonzero color electric field. Thus for the semi-QGP, one 
does not expect a nonzero color field in vacuum, and such terms should not appear. We 
admit that at present, we do not have a fully self-consistent theory of the semi-QGP, which 
would allow us to demonstrate this. 

We turn to subleading terms in the gluon self-energy, ~ T^. In this, it helps greatly to 
recognize that the integral must be even in x = cos^^, Eq. (110). The measure is even in x, 
as are terms ~ k^k\ which produce contributions ~ 1 or ~ x"^. 

There are several ways that corrections ~ can arise. The first is a numerator propor- 
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tional to and p' : 
1 



327r2 



dk k {n{k — iQi) — n{k + iQi) + n{k — iQ2) — n{k + iQ2)) 

(123) 



4vr V / \ip^lp — f)-k ipl'^/p + p-k 

This vanishes, because the integrand is odd in k. The second is expanding 1/Ep_k which 
arises in the measure of the integral, as the residue of the propagator, 

1 1 p ■ k , , 

+ V + (124) 



Ep^k k k'^ 

Again, since p-k = px, this is odd in x, and so vanishes. 

The third is by expanding Ep^k in the energy denominators: 

2 



1 1 p^-[p-k 

T—^ (125) 



tpl^ ± {Ek - Ep.u) ip}^±p- k 2k (tpl^ ±p-ky 

The numerator of the second term on the right hand side is — x^), which is even in x. 
Because of the =F sign in front of the second term, though, this is in all odd in x, and so 
vanishes. 

Thus, the only way that corrections ~ arise is by expanding Ep_k in the statistical 
distribution functions 

n{Ep_k - iQ) ^ n{k - iQ) + ^ n{k - iQ) (1 + n{k - iQ)) + ... (126) 

This is exactly the same sort of expression as at Q = 0. Now this term is odd in x, but note 
that —n{Ep_k + iQi,2) enters inX2, Eq. (116), and +n{Ep_k — iQ 1,2) enters in X3, Eq. (117). 
Thus the result is even in x, as it must be if not to vanish. Explicitly, the terms ~ in 
J"^^ are a product of an integral over k, 

1 

dk k'^ {n{k - iQi){l + n{k - iQi)) + n{k + iQi) (1 + n{k + iQi)) + {Qi ^ Q2)) ■ 

(127) 



167r2T 



and an integral over the angular variables, 

h^h.l r) ■ k \ 

(128) 



dn / k'k^ p- k 
An \ p.K 
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The angular integral in Eq. (128) is identical to that for ordinary hard thermal loops. 
The momentum integral is also a minor modification. Consider 

- / dk P {n{k - iQ){l + n{k - iQ)) + n{k + iQ){l + n{k + iQ)) , (129) 







which we can rewrite as 

-2—/ dk k"^ {n{k - iQ) - n{k + iQ)) . (130) 

This integral arose previously in Eq. (102), and involves the function Ao{Q)^ which is the 
first derivative of the potential V{Q), Eq. (95). Since in Eq. (130) we take a derivative of 
this function with respect to Q, however, the momentum integral in Eq. (129) involves not 
the first derivative of V{Q), but the second, through the function A{Q), Eq. (96): 

AQ) ■ (131) 
In summary, the terms ~ in J"^^{P, Qi, Q2) are those of Eq. (120); the terms ~ are 

T2 ^ ,^ , S'^ f dVL k^y p- k 



(.4(Q0+.4(Q.))^- + y-^J . (132) 

For completeness, we have added the tadpole terms, ~ 6"^^ /Q. This expression can be rewrit- 
ten as 



- {Am + Am) y - j ^ . (133) 



It is direct to show that these are the only terms ~ or ~ T^. The one concern is 
terms where the numerator is ~ k'^p' or ~ p'^k^: while the single power of brings in a 
suppression by a soft momenta, ~ p/T, this times a term ~ T^/p could produce a result 
~ T^. However, a single power of fc* or y is manifestly odd in the angular variable x, 
Eq. (110). In contrast, the terms which produce contributions ~ T^, Eq. (109), are even in 
X. Thus these possible terms ~ vanish when integrated over x. Therefore, in the hard 
thermal loop approximation, terms proportional to and pP in Eq. (Ill) may simply be 
dropped, 

^ ^ ^ HTL /• d^K k' y , , 

For the gluon self-energy, the general integral required is 
J^^P'^ Q„ Q,) =\ I^A- ( Tl^TV^ + ^ 



4 J (27r)4 V \{K + Qi)^ (K + Q^)^ 

{2K-P + Q1- Qa)^ {2K-P + Q1- Q, 



+ 



{K + Q,)\P-K + Q2)^ 
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(135) 



where = u^Q, m'^ = 1 if /i = 0, and zero otherwise. The first term on the left hand side, 
~ 5^^ , is added in anticipation of the integrals which arise for hard thermal loops, and is 
obviously independent of the external momentum, P. 

Momentum dependence arises from the second term on the left hand side. Its numerator 
involves 2K — P + Q1 — Q2 = "^K^ — P^^] under a shift of the loop momentum, K —>■ P — K, 
this becomes —(2K'^ — P^^). Comparing with Eq. (112) shows that Eq. (135) is symmetric 
under interchange of Qi and Q2. 

When Q = 0, the computation of hard thermal loops is simplified by keeping only powers 
of the hard loop momentum in the numerator, and dropping those of the soft external 
momenta, P'^ [23, 24, 25]. When Q ^ 0, the analogous approximation is to keep powers of 
or K^, and drop powers of P^^, as we saw in J'^K It is essential to remember that while 
Qi and Q2 are each hard, that the real energy P^^ is soft. Sec. IIIC. Thus we can write the 
second term on the left hand side of Eq. (135) as 

d^K {K^Y {K^y HTL r d^K {K^Y {k^y 

(2^ {K^Y{P^^ - K^Y J (2^ {K^Y{P^^ - K^Y ' 
This simplification will help greatly in computing the gluon self-energy in Sees. VD and 
VE. 

Computations similar to those above show that this integral equals 

J'''{P'\Qi.Q2) ^ (A(gi) + A(g2))5r^'^(P'')+^ {A{Q,) + A{Q2Y5n^''{p'^). 

o 24 

(137) 

where the momentum dependence enters through the functions Sr^'^{P^'^), Eq. (121), and 

*J7-(P) = (^-uV -^p,j'^ . (138) 

Again, K = {i, k), so K"^ = 0, and the angular integral is over all directions of k. 

Previously we computed SP'^^ and 611^^ . Given this result, it is immediate to show that 
6r^^ and SP^^ are correct. The only effort necessary is to establish the correctness of 5i7*° 
and 5n^^. However, the derivation of 511^^ above shows that the terms ~ arise in 
precisely the same way when Q 7^ as for Q = 0, entirely from the change in the statistical 
distribution functions, Eq. (126). Thus we can be certain that the only change is that of 
the thermal gluon mass in a background field. That is, the momentum dependence, which 
enters through 611'^'^ (P), is unchanged. For this, the previous computation suffices. 
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The function 6n^^"{P) is the same function as arises in the hard thermal loops when 
(5 = 0. The function 5r^^^{P) enters only when Q 7^ 0, since it is multiplied by Aq, and 
^o(O) = 0. From Eqs. (121) and (138), the two functions are simply related to one another, 

^r^^lP) = — (Sn'^'iP) + u^v") . (139) 

The function 6n^'^{P) is transverse, 

sn>"'{P) = . (140) 

Thus the new hard thermal loop is not, 

P" Sr^'^iP) = iu" . (141) 

As discussed following Eq. (122), the new hard thermal loop is the modification of the 
propagator in a background field. Thus it is not necessary for the term to be transverse. We 
shall also see this in the next section, where we compute the one point function for a gluon. 



C. One point gluon function 

With these results in hand, the diagrams can be computed directly. In this subsection we 
begin with the one point function for a gluon. The integral which enters is that of Eq. (100), 
and generates the function Aq, Eq. (102). As can be seen from Eq. (101), this function 
vanishes when Q = 0, as then the integral is odd in k^. 

Consider a gluon loop tied onto a gluon line. By an argument similar to that which lead 
to Eqs. (50) and (51), if the gluon propagator in the loop is 'P'^'^'^f / (^K'^^Y , Eq. (39), we can 
replace this by 5^^ 5'''^ / {K'^'^Y . Since the external gluon has zero momentum, the three gluon 
vertex involves one momentum, which we choose as i^^'^; from Eq. (45), the three gluon 
vertex is 

- {-K^ 5^, + 2K';f 5,, - Kf Sx,) . (142) 

This is transverse in K^"^ and K'^f, so the gauge dependent term in the gluon propagator, 
~ (^ — 1)K^'^K^'^, Eq. (39), drops out. The term in the gluon propagator ~ 5^*^ acting upon 
Eq. (142) gives 2{d — 1) in d = 4 spacetime dimensions. Including an overall 1/2 for the 
gluon loop, 
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Here = (1,0). The contribution from a ghost loop is similar, with a coefficient of — 1 
instead of +3 in Eq. (143); thus the sum is +2, which reflects for the two transverse degrees 
of freedom of a gluon. 

The contribution of A^f flavors of massless quarks is 

(-)'9N, (r^)^^ J j2^^''Z[f-c = Y«'^ Mcc MQ') ■ (144) 

Using Eqs. (7) and (24), the sum of the gluon, ghost, and quark contributions is 




NfAoiQ") ■ (145) 



In obtaining this, we have used the fact that is odd in Q — > —Q- This is clear from its 
definition in Eq. (101); its form in Eq. (102) is only valid for < g < 1. 

For a Z(A^) interface, this current is part of the equation of motion for the gluon, corrected 
to one-loop order. There it is natural: the interface arises from a balance of the Lagrangian 
at tree level, and a potential for the Q's at one-loop order. Equation (145) is precisely the 
derivative of the one-loop potential. 

For the semi-QGP, a term must be added to the Lagrangian to cancel the contribution of 
this current. This is natural; if such a term is not added, the minimum would be at Q = 0, 
or equivalent points; i.e., the usual Z(A^) minima. 



D. Quark contribution to the gluon self-energy 

With the previous examples in hand, the computation of the gluon self-energy is mainly 
a matter of putting things together. Even so, because the Q's are nontrivial, what is of 
interest is to see how new color structures arise at one-loop order from quantum corrections 
to the propagators at tree level. This is unlike the case of the quark self-energy, Eq. (87), 
where the only change from Q = was the value of the thermal quark mass. 

We start with the contribution of the quark loop. For A^f flavors, this is 

- (^;r).k(P"^) = (-) {^9rN, f^.tr ^ {e% -L- 7^ it'')fe—~^— ■ 

(146) 

Keeping only the terms in the numerators, the hard thermal loop in this contribution 
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FIG. 9: One-loop diagram for the quark loop in the gluon self-energy. 



IS 



ej=l 



N 



(147) 



Because Q*^ appears in the function J''^", what enters is not simply the trace of two projection 
operators, Eq. (8). The sum over the color index e and / is easy, and gives 



T-TTT 



1 



N 



1 ^ ^ ^ 



e=l 



This is illustrated in Fig. 9. There is the usual planar diagram, plus three contributions 
from diagrams in which one or both of the gluon indices are traced. Notice that as usual, 
the color structure is far more clear from the diagram, than from the detailed expression in 
Eq. (148). 

The last three terms, ~ are like the photon propagator for QED, where the 

fermions propagate in a background Q field. In all, the quark contribution is traceless, 

N N 

J2 {n;r%dP) = J2^nfrMP) = . (149) 



a=l 



c=l 
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which is necessary for self-consistency. Note that the quark contribution is not simply 
proportional to a projector operator, ~ -pa^M^ Instead, because the Q's can be unequal, at 
one-loop order the color structure is more complicated, Eq. (148). 

E. Ghost and gluon contributions to the gluon self-energy 

To compute the contributions of gluons and ghosts to the hard thermal loop in the gluon 
self-energy, we first compute in background Feynman gauge, ^ = 1 in Eq. (39). We then 
show that the results are independent of ^, in complete analogy to when (5 = 0. 

The simplest contribution is the tadpole diagram, which is independent of the external 
momentum. This involves the four gluon vertex of Eq. (49). For the gluon in the loop, we 
can use the simplification of Eq. (51) to take the gluon propagator as a Kronecker delta. 
The result is 



The coefficient 6 2{d — 1) m. d spacetime dimensions; the 1/2 is for a bosonic loop. 

The gluon and ghost loops are, generally, involved. However, we can use the simplification 
of Eqs. (136): write the momenta in terms of the external momentum, P""^, and a loop 
momentum, which we can define as K^'^. Then although the Q's are hard to begin with, we 
can consistently treat P"** as soft, and K^^ as hard. This allows us to drop powers of P*^* 
uniformly. (This is only valid in Feynman and Coulomb gauges [23, 24, 25]). 

With this approximation it is then easy to read off the hard thermal loops in the gluon 
self-energy. The ghost loop is given by taking K^^ at each vertex, and so is 



For the contribution to the gluon self-energy from the diagram with two three gluon vertices, 
in each vertex we can neglect the external momentum, as in Eq. (142). This gives 




(150) 




(151) 




[2^ {Kf^YiP'^'^ - Kf- f 



(152) 



The coefficient 10 2{2d — 3) in d spacetime dimensions. 
In all, the sum of Eqs. (150), (151), and (152) is 



(i7^t'^')gi(^"') 



HTL 



,2 r{ab,ef,gh) £{cdje,hg) j^iu/pab Qfe Qhg 



(153) 



5V ;V 



41 



d- 
c- 





.a 
-b 




j^fe _ pab 



FIG. 10: One-loop diagram for the gluon loop in the gluon self-energy. 



From Eq. (135) 



FTTT 

ab,cd\ ^ ( pab\ '^^^ 



(154) 



Like the quark self-energy, Eq. (149), this is traceless, as it must be, to represent a matrix 
in SU(A^). 

The color structure is illustrated in the diagram of Fig. 10. There is the planar diagram, 
minus a diagram in which the indices are summed over. 

In all, the hard thermal loop in the gluon self-energy is the sum of Eqs. (148) and (154). 
Each term is a product of a color dependent factor, times a function of the soft momentum: 

HTL 



The thermal gluon mass in a background field is 



(155) 



ad jrbc 



6 



5"% 



N 



-2 S^'d' 



ab jrcd 



e=l 

Nf ( 



In zero field, ^(0) = 1, AiQ'') = A{nT) = -1/2, and 



ab,cd 



N + 



2 y 3 

The new hard thermal loop in the gluon self-energy involves the color matrix: 



(157) 



2n-^3 xo.d xbc 



N 



{MQn + AoiQ'^')) - N, ( A(Q°) + A(-Q') 



e=l 



(158) 



The terms proportional to 6°'^6'^'^ vanish because Ao{Q) is odd in Q, Eq. (102). As necessary, 
the matrix }C°'^'^'^{Q) is traceless. 
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Comparing Eq. (158) with Eq. (145), we find 

This is natural; the background field induces a color current, which couples to the gluon. 
The self-energy obeys the Ward-Takahashi identity in the background field: 

These expressions were computed in Feynman gauge, but the results are independent of 
the gauge fixing parameter, ^, Eq. (29). Except for the gluon self-energy, where the tadpole 
diagram enters, the hard thermal loops in any gluon amplitude only involve three gluon 
vertices. This vertex satisfies an Abelian-type Ward identity, Eq. (47). This identity can 
be used to show that all hard thermal loops are independent of ^, up to possibly gauge 
dependence in terms which are independent of the static momentum. However, it is known 
that at one-loop order that the potential in a background Q field is independent of C, [6, 7, 8]. 
Thus, the hard thermal loops are independent of ^ when Q ^ 0, as for Q = 0. 

VI. CONCLUSIONS 

In this paper we developed techniques to analyze the real time response functions for a 
't Hooft loop, or Z(A^) interface. By introducing the double line notation in Sec. II, we are 
able to analyze a much more general problem, as is appropriate for the semi-QGP phase of 
a gauge theory [32]. 

While the final expressions which we obtain appear involved, in fact the physics for Q ^ 
is very similar to that for Q = 0. For the quark self-energy, Eq. (87), the hard thermal loop 
is a thermal quark mass times a function of momentum. The function of momentum, 5S(P) 
in Eq. (88), is unchanged from Q = 0. What does change with Q is the thermal quark mass, 
Eq. (89). It is most natural that in the presence of a nonzero background field, that the 
curvature about the minimum changes with the background field. 

For the gluon self-energy, Eq. (155), there is a piece very similar to that for Q = 0. There 
is the same function of momentum, 511^'^ {P) in Eq. (138), as in zero field. This function is 
multiplied by a thermal gluon mass, Eq. (156), which is of course Q-dependent. 

The surprise is that there is a new function in the gluon self-energy, Sr^'^{P) in Eq. (121). 
(Note, however, that this function is linearly related to 611^'^ {P), Eq. (139)). The usual 
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hard thermal loop in the gluon self-energy is SU'^" ~ g'^T'^, and is smaller than the new hard 
thermal loop, ST^'^ ~ g'^T^/p. The general principle of hard thermal loops, however, is that 
for soft external momentum, they are as large as the terms at tree level. This remains valid, 
since for a non-Abelian gauge field, the propagator in a background field has a new term, 
Eq. (30). 

We also note that surely our entire derivation would be much simpler if we had used 
kinetic theory in the presence of a nonzero background field. However, we preferred to use 
an ordinary perturbative analysis, since computing at Q 7^ is an unfamiliar exercise. We 
do expect that the derivation of a complete action for all hard thermal loops for nonzero Q 
would be much simpler with kinetic theory, as it is when Q = 0. 
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APPENDIX: Z{N) AND U(l) INTERFACES 

In this appendix we use this opportunity to make some comments about Z(A^) [3, 4, 5, 6, 7] 
and U(l) [29] interfaces in both Abelian and non-Abelian gauge theories. 

As discussed in the Introduction, Sec. I, to define an interface we pick out one of the three 
spatial directions, say that in the z direction, and consider it separately from the two other 
spatial directions, and from that for imaginary time, r. It is also necessary to assume that 
the length in the z direction, L, is much larger than that for the two transverse dimensions, 
of size Lt [6] . Lastly, we assume that both L and Lt are much larger than any physical mass 
scale, such as the inverse Debye mass, ~ l/{gT). 
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For a pure gauge theory, without dynamical quarks the gauge group is SU(A^)/Z(A^). 
Thus we can require that the gauge field at z = L is a gauge transformation of that at 
z = 0. We take this gauge transformation to be a constant element of Z(A^) [3]: 

(^ = exp — 1^. (A.l) 



N 

Since this VL commutes with all group elements, the gauge field at 2; = L is identical to 
that at z = 0. In going from z = to z = L, though, one winds, in a topologically 
nontrivial manner, in the configuration space of gauge potentials [44]. One can show that 
these boundary conditions are equivalent to inserting a 't Hooft loop, at z = L/2, at the 
boundary of the two transverse dimensions [3, 7]. 

There are k distinct transformations possible, where k = 1,2, . . . {N — 1). This generates 
the cyclic group, 'L{N), where k = N is equivalent to the identity. At tree level, this 
transformation is implemented by a gauge transformation which is linear in z, 

^{z) = e^^{2mk (y2t^^^ |) . (A.2) 

The factor of \/2 is because of the normalization for t^^ in (10). This gauge trans- 
formation is generated by a constant vector potential in the imaginary time direction, 
Aq = {2'KikT / g){\/2t^^)z / L. This is a linear ansatz, and is only valid for the tree level ac- 
tion. Loop corrections generate a potential for Aq, and turn the true solution into a domain 
wall, of width ~ l/{gT) [6]. At present, however, we need not trouble ourselves with such 
details. All that matters is that this can be modified to represent a solution of the quantum 
equations of motion. 

Consider the sum of the following two diagonal generators: 

A3 = 72 {{N - 1) + 1^^) = ^2{N-l){N-2) A^^i , (A.3) 

Eq. (11). 

Consider a gauge transformation generated by this transformation, Q{z) = 
exp(27ri A3 z/L)!^. This is a "U(l)" transformation, in that it is the same on both ends 
of the box, VL{L) = Q{0) = In [29]. This is evident from (A.3), where we see that it is a 
combination of two Z(A^) transformations: one of strength — 1, along the t^-i-^-i direc- 
tion in group space, plus one of unit strength, along the direction. The total strength 
is then A^, which for Z(A^) is equivalent to the identity; i.e., to no winding. 
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In [29], U(l) interfaces were suggested as a way to perform the semiclassical matching 
between an effective and an original theory. For the pure gauge theory, the U(l) interfaces 
are not necessary: we can use the k possible Z(A^) interfaces to perform the matching. We 
remark, rather trivially, that if the Z(A^) interfaces do match, then so will the U(l) interfaces, 
since by the above, they are just a combination of two Z(A^) interfaces in different directions 
in group space, one along t^^^'^^^, and the other along t^'^ . 

The real use of a U(l) interface, however, is for the theory with dynamical quarks, where 
the Z(iV) invariance is broken by the presence of quarks. Before considering this case, 
consider a simpler example: a Higgs model, where a charged scalar field, 0, acquires a 
vacuum expectation value, (0). Then in going from one end of the box, to the other, one 
can define a winding number by the number of times that the phase of (0) winds around 
2tt. The topology is elementary, just 7ri(U(l)) ~ Z. 

Now consider an Abelian gauge theory at nonzero temperature, such as QED, where the 
gauge group is unbroken. Instead of the expectation value of a Higgs field, we consider 
windings of the thermal Wilson line. The thermal Wilson line is = +1 at both ends of 
the box, but it can wind nontrivially as it goes along the box. The topology remains 
7ri(U(l)) Z. Thus there are thermal interfaces in QED. This was first pointed out by 
Smilga [15], who computed their properties in the Schwinger model. 

Thermal U(l) interfaces are of interest in the electroweak theory. They are unlike stan- 
dard domain walls, in that the interface tension is strongly dependent upon the temperature. 
In the high temperature phase, where the Higgs field is unbroken, walls lie in the U(l) for 
hypercharge. As the system cools through the electroweak phase transition, they then rotate 
into the U(l) for electromagnetism. Since the potential is generated by quantum effects, it 
vanishes exponentially at low temperatures, exp(— me/T), where rrie is the electron mass. 
Thus while U(l) domain walls dominate the stress energy tensor while they exist, unlike 
standard domain walls, they naturally vanish at low temperature. Whether their presence, 
during the period in which they dominate the stress energy tensor, can lead to characteristic 
cosmological signals is an interesting question. 

We return to a SU(A^) gauge theory. For each of the diagonal directions, we can 
define a U(l) winding number. For example, A3 represents winding once in the first N — 2 
directions, — (A^ — 2) times in the (A^ — 1)*^ direction, and no winding in the last, A^*'^ 
direction. 
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Semiclassically, one expands about the configuration along the A3 direction, in both the 
original and effective theories. The configuration is a local minimum of each action, and has 
no instabilities under small fluctuations. Thus such a U(l) interface can be used to match 
the parameters of the effective theory to the original theory, as suggested in Ref. [29]. 

To define such a configuration nonperturbatively, such as on the lattice, it is necessary 
to fix the freedom to perform global gauge rotations. Only after fixing the freedom to 
change overall gauge rotations can one define U(l) winding numbers for each of the diagonal 
directions [29]. 
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